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1.  INTRODUCTION 


Currently,  the  most  widely  used  design  for  kinetic  energy,  antitank  applications  is  the 
gun-launclu'd,  fin-stabilized,  long-rod  projectile.  The  cross-sectional  diameter  of  the  rod  is 
much  smaller  than  the  diameter  of  the  gun  bore.  Fins  span  the  area  bet  .vi'en  the  rod  and  the 
gun  tube.  Therefore,  a  sabot  (or  carrier)  is  required  to  provide  obturation  for  the  i)rojr‘ctile 
and  to  minimize  its  in-bore  balloting.  Once  free  of  tlie  gun  tube,  tlie  sabot  must  be  discarded 
in  order  to  permit  unconstrained,  low-drag  flight  to  the  target.  The  sabot  is  divided  into 
three  or  four  components  along  axial  planes.  For  smooth  bor<'  gun  tubes,  these  comi)onents 
separate  from  the  projectile  under  the  action  of  elastic  and  aerodynamic  loaib.  Figures  1 
and  2  show  a  photograidi  and  sharlowgraph  of  typical  sabot  discard  during  free  llight. 

It  has  been  demonstrated  (Schmidt  el. al.  1978)  that  aeiodynamic  intiuference  geiK'rateil 
by  the  sabot  components  can  be  a  significant  source  of  projectih*  launch  disturbanct'  h'ading 
to  unacceptable  lo.s.s  of  accuracy  at  the  target.  Perturbations  to  the  projectile's  trajt'ctory  are 
magnified  by  geonu'tric  asymmetry  in  the  discard  pattern  and  Iry  extended  periods  during 
launch  when  the  sabot  components  are  in  close  proximity  to  the  projectile.  .A  detailed 
understanding  of  the  three-dimensional  shock/boundary- layer  intcrf(>rence  flowfleld  between 
the  .saf)ot  and  the  projectile  (see  Figure  2)  is  not  available. 

.An  extensive  experimental  program  to  inve  ligate  tlie  aerodynamics  of  sabot  discard  has 
been  conduct<*d  (Schmidt  1981).  During  these  tests,  a  proji'ctilc  and  three  sabot  components 
were  sting-mounted  in  the  N.ASA  Langley  Unitary  Plan  wind  tunnel  facility  1  x  I  ft  test  sec¬ 
tion.  The  model  configuration  included  a  stationary  cone-cylinder  projectile  (without  fins) 
at  zero  angle-of-attack  and  three  120°  inclnded-angle  sabot  components  located  symmetri¬ 
cally  around  the  projectile.  Figure  3  shows  a  schematic  (cross-section)  of  the  wind  tunnel 
model  (one  sabot  shown).  The  cylinder  section  of  the  projectile  was  50.8mm  in  diaiiu'ti'r;  the 
projectile  had  a  length-to-diameter  ratio  of  10.5  and  a  30°  included-angle  conical  nose.  I'ifly 
static  pressure  taps  were  positioned  on  the  surface  between  the  120°  planes  of  symmetry, 
with  four  taps  on  the  conical  section.  The  sabot  had  cylindrical  inner  and  outer  surfaces 
of  radii  25. 1  and  76.2mm,  respectively,  with  the  a  leading  edge  chamfer  of  10".  Fifty  static 
pressure  taps  were  located  on  the  inner  and  outer  surfaces.  The  test  Mach  number  and 
Reynolds  number  were  4.5  and  6.G  million  per  meter,  respectively.  A  typical  flight  Reynolds 
number  of  89  million  per  meter  could  not  be  reproduced  in  the  tunnel.  Test  res\dt.s  showed 
regions  of  shock/boundary- layer  interaction,  separated  flow  and  other  viscous  phenomena 
which  are  sensitive  to  the  Reynolds  number. 
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Initial  analytical  work  for  steady  state  sabot  discard  aerodynamics  udied  on  Uic  Newto¬ 
nian  How  approximation  and  empirical  aerodynamic  interaction  analyses,  for  cxampb',  the 
.\\('()  cocie  (Ciimi  ct.al.  1977.  Seigelman  et.al.  IH^d).  (\inside,ation  was  limiicd  to  a  gen¬ 
ii  al  sabot  coniignrat  ion,  bound  radially  by  two  cylindrical  siirfai  es  and  axially  l)\  two  conical 
siiilaces.  riiese  assnm|)tions  make  discard  coin()ntaf ions  tractable  and  in  some  cases  rej) 
resent  ac< urate  apirro.ximations.  However,  it  is  apparent  that  the  multiple  shoe  k/exi>ansion 
interaction  llovvtii'ld  betweeti  the  projectile  and  sabot  petals  is  ati  essimtial  |)art  of  the  anal¬ 
ysis.  rite  initial  version  of  the  .NVl'O  code  (Criini  et.al.  1977)  evaluated  the  aerodvnamic 
loadings  on  the  sabot  .segments  using  Newtonian  theory  and  a  snbsotiic/suirersonic  inlet 
model;  jtressure  forces  on  each  surface  of  the  segments,  indnding  sabirt  siries,  were  obtained 
separately  and  summed  to  provide  results  for  total  force  and  moments  (excluding  shear  st  ress 
( um|)onent  s).  llu'  code  assumed  that  the  aerodynamic  coedicient  s  lor  the  projectile  were 
known.  .Although  the  sabot  separation  proci'ss  is  initially  dominated  liy  aerodynatmc  interac- 
I  ion.  the  code  assumed  otie-dimensional  How  between  the  botlies.  Recent  versions  (Seigelman 
et  .al.  19S,{)  include  an  integrated  How  elcnnent  approach  utilixing  local  shock/c'.xpansion  pro- 
cedures  bitseni  on  saliot  surface  pressures  measured  during  wind  tunned  tests  (Schmid'  19S1). 

I  hese  test  data  are  used  to  determine  pressure'  ic'vels  on  certaiti  sabot  locations  with  line'ar 
\atiatiotis  assumed  betwe'on  these  points.  As  a  reesult,  the  cerde  include'S  the  effects  of  pres¬ 
sure'  pulses  ott  tile  bodies  caused  by  impinging  and  re'flecting  shock  waves.  'Id  calculat.*  the 
initial  sabot  lift-off  aerodytiamic  force's  and  then  again  when  the  sabot  pc'tals  are'  not  in  close 
proximity  to  the  projectile,  Ne'wtoniari  flow  the.'ory  is  used.  In  some'  cases,  howe*ve>r,  these 
code  improvements  produced  overestimates  of  the  discard  proee'ss  in  contrast  to  initial  code 
jiredict  ions. 

1  his  report  describes  computational  fluid  dynamics  (CFD)  .solutions  applied  to  the  three- 
dinu'nsional  (3D)  Navier-Stokes  equations  for  symmetric  sabot  discard.  During  symmetric 
discard  multiple  sabot  components  are' assumed  to  follow  identical  trajectories  away  from  the 
proje'ctile,  and  the  projectile  is  assumeel  to  be  at  zero  angle-of-attack.  As  shown  in  Figure 
1.  the'  eeunputational  domain  can  therefeue  be  limited  to  a  smaller  portion  of  the  entire 
flowfie-ld  around  the  configuration;  this  reduems  computational  grid  si/e,  computer  memory, 
and  computer  run  time.  For  three  sabot  components  this  domain  spans  a  bO"  sector  from 
sabot  midplane  to  symmetry  plane  betw'een  neighboring  sabot  components.  For  asymmetric 
discard  the  computational  domain  would  be  greatly  expanded  (i.e.  a  full  360"  se'ctor)  with 
a  corre'sponding  increase  in  computer  requirements.  The  portion  of  the  launch  cycle  that 
involves  strong  aerodynamic  interference  between  the  projectile  and  the  sabot  components 
is  examined.  Thus,  simulations  are  performed  for  small  vertical  separation  of  the  sabot  from 
the  projectile  surface.  Ay/D  <  1  (D  =  projectile  rod  diameter  =  1  cal.  or  50.8mm  in 
l  igiiri'  3)  and  sabot  angle  of  attack  o  f"  10".  Previous  work  de.scribed  code  validation  with 
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wind  tunnel  results  (Nusca  Aug.  1990,  Oct.  1990).  A  four-stage  sabot  discard  sequence  was 
numerically  simulated  for  the  wind  tunnel  model  configuration  (Nusca  Apr.  JuI.  1991). 

In  this  report  these  simulations  have  been  extended  to  ten  stages  with  resultant  aerodynamic 
forces  and  moments  computed  from  the  flowfield.  The  symmetric  sabet  discard  trajectory  is 
then  simulated  and  compared  to  results  obtained  using  the  AVCO  code.  This  quasi-steady, 
programmed  simulation  ignores  the  flow  time  dynamics  and  does  not  link  the  aerodynamic 
forces  to  the  sabot  motion.  However,  such  a  simulation  serves  as  a  prelude  to  computations 
that  utilize  coupling  of  unsteady  aerodynamics  and  rigid-body  motion.  The  flowfield  for  a 
M865  projectile/sabot  is  also  simulated. 

2.  COMPUTATIONAL  APPROACH 

('I'D  can  be  used  to  simulate  the  compressible  flowfield  around  aerodynamic  bodies  bv 
solving  tlie  3D  Reynolds-averaged  Navier-Stokcs  (RANS)  equations.  The  US.A-PGfi  code  was 
dev('loi)ed  l)y  ('liaki avarthy  (1985,  1988).  The  R.-XNS  equations  are  written  using  a  perfect 
gas  assumption.  Both  laminar  and  turbulent  flows  can  be  investigated  thus,  a  turbulence 
model  (Baldwin  1978)  is  required  for  closure.  In  addition,  backflow  regions  can  be  present 
thus,  a  backflow  turbulence  model  (Goldberg  1986)  is  included.  The  equations  are  trans¬ 
formed  into  the  conservation  law  form  and  discretized  using  finite- volume'  api)re)ximations. 
The  I'S.-X-PGfi  code  uses  a  class  of  numerical  algorithms  ternu'd  total  variational  diminishing 
(TVD).  I'he  resulting  set  of  ('quations  is  solved  using  an  implicit,  factored,  time-stepping 
algorithm.  The  solution  takes  place  on  a  computational  grid  that  is  generated  around  the 
configuration  in  zones;  zonal  boundaries  are  transparent  to  the  flowfield. 


2.1  Equations  of  Motion.  The  RA.N.S  equations  for  3D  flow  an'  written  in  the'  fol¬ 
lowing  conservation  form.  I’lie  dependent  variables  u,  v,  w,  and  c  are  mass-averaged. 


dW 

ur 


OF  dc  di{ 


(1) 


ir  =  (p  f)u  pv  pir  pc) 


(  pii  ^ 

pu^  - 
par  -  Tj.y 
piiir  —  Tj - 

\  put  +  (),  -  rj,,iL  -  t,,w  / 
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Arrays  (!  and  H  are  similar  in  form  to  array  F  (see  Nusca  Oct.  1990).  Normal  stress 
(rr),  shear  stress  (r),  heat  transfer  [q]  and  energy  (e)  arc  defined  elsewhere  (N'nsca  Oct. 
1!>90).  d'he  laminar  and  eddy  viscosities,  ^  and  //<,  are  implicitly  divided  by  tlie  reference 
Keynolds  number  (He).  The  flow  medium  is  assumed  to  be  a  perfect  gas  satisfying  the 
eipiation  of  state  p  =  p^T.  A  power  law  (Mazor  et.al.  1985)  is  u.sed  to  relate  molecular 
viscosity,  p,  to  temperature.  'I'he  laminar  and  turbulent  Frandtl  numbers.  Pr  and  Pr,,  are 
assumed  to  be  constant  with  values  of  0.72  and  0.9,  respectively.  The  ratio  of  specific  heats, 

7,  is  also  assumed  constant  Assuming  a  time-invariant  grid  and  using  the  transformation  of 
coordinates  imidied  by  r  =  <,  ^  =  q  =  r/(x,  (/,;■),  and  (  -  Pcpiation  1  ran 

be  recast  into  the  conservation  form  where  q,  and  (,'  are  the  new  indepeiuh'iit  variables, 
and  .1-^,  Xr,,  X(,  y^,  y,,,  and  Z(  are  the  nine  transformation  coefficients  obtained 

numerically  from  the  mapping  procedure.  Transformed  time  is  represented  by  r. 

-TT  f  T  [{yt)^  -  +  f'V.v  “  f^/y]  =  o  (2) 

OT  .'\rca 

I  he  “Area"  in  Equation  2  denotes  the  area  of  the  finite  volume  cell  being  considered  at  the 
time  of  discretization  of  the  equations. 

The  shock/boundary-layer  interference  flowfield  between  projectile  and  sabots  can  in¬ 
clude  regions  of  recirculating  flow.  To  improve  the  predictive  capability  of  separated  flows 
using  HANS  codes  a  new  turbulence  model  has  been  recently  developed  by  Goldbi-rg  (1986). 
The  new  model  is  based  on  experimental  observations  of  detached  flows  and  allows  turbulence 
due  to  local  shear  effects  to  be  taken  into  account  in  addition  to  wall-shear  contributions. 

'I  he  velocity  scale  function,  which  is  normally  yu?,  is  modified  as  [y  —  yf.)^'  (for  y  >  y^).  Here, 
u.’  is  the  magnitude  of  the  local  vorticity  and  y^  is  the  location  away  from  the  wall  where  the 
vorticitv  first  diminishes  to  a  small  fraction  (<  5%)  of  the  local  maximum  magnitude.  hTom 
this  location  onward  the  length  scale  is  given  by  y,„ax  —  y^.  The  model  prescribes  turbulence 
kinetic  energy  and  dissipation  analytically  within  backflows.  A  formula  for  the  eddy  viscos¬ 
ity  (p,)  within  backflows  is  derived  and  used  for  the  HANS  equations  when  calculations  are 
done  inside  separation  bubbles.  Outside  of  them,  another  turbulence  model  (Baldwin  et.al. 
1978)  supplies  the  values  of  eddy  viscosity. 

2.2  Computational  Algorithm.  The  spatial  discretization  techni<iue  for  the  e()ua- 
tions  of  motion  must  successfully  capture  the  complex  physics  of  interacting  projectile/sabot 
flowfields.  '[lie  'FVD  formulation  for  the  convective  terms  along  with  a  .special  treatment  of 
lh('  dissijiat  i  v('  terms  (Eaiuation  1)  provides  an  appropriate  simulation.  In  recent  years,  d'VI) 
formulat  ions  have  been  constructed  for  shock  ca|)turing  finite-difference  nu't  hofls  ( ( 'hakravart  hv 
et.al.  1985,1988).  Near  large  gradients  in  the  solution  (extn'ma)  lA''!)  algorithms  automat¬ 
ically  reduce  to  first-order  accurate  discretizations  locall\'  while  away  from  extriuna  they 


A 


can  be  constructed  to  be  of  higher-order  accuracy.  This  local  effect  restricts  the  niaximuni 
global  accuracy  possible  for  TVD  algorithms  to  third  order  for  steady-state  solutions.  TVl) 
methods  manifest  many  properties  desirable  in  numerical  solution  procedures.  By  design 
they  avoid  numerical  oscillations  and  ‘‘expansion  shocks”  while  at  the  same  time  being  of 
higher-order  accuracy.  TVD  formulations  are  also  based  on  the  principle  of  discrete  or 
nunu'rical  conservation  which  is  the  numerical  analog  of  physical  conservation  of  mass,  mo¬ 
mentum,  and  energy.  Thus,  TVD  algorithms  can  “capture”  flowTield  discontinuities  (e.g. 
shock  waves)  with  high  resolution.  At  a  fundamental  level  they  are  based  on  upwind  al- 
goritlmis;  therefore,  tlu'y  chjsely  simulate  the  signal  propagation  properties  of  hyi>erbolic 
etiuations.  Algorithms  based  on  the  TVD  formulation  arc  completely  delined.  In  contrast, 
algorithms  based  solely  on  central  differences  involve  global  dissipation  terms  for  stability 
and  have  one  or  more  coefficients  that  must  !)e  judiciously  chosen  to  achieve  desirable  re¬ 
sults.  .Any  conventional  lime  discretization  nu'thod  suitable  for  the  .\avier-Stokes  ('(juations 
can  b(>  used  togetlu'r  with  this  space  discretization  nu'thodology :  for  examph',  ai)|)roximati‘ 
factorization  and  r<-laxation  t(‘chni<)ues. 

2.3  Computational  Grid.  Numerical  simulation  of  the  interacting  flowfield  about 
projectile/sabot  combinations  is  complicated  by  the  non-axisymmetric,  multijile-wall  geom¬ 
etry.  d  lx.-  computational  domain  is  divided  into  zones  of  simjile  geometric  shape.  In  each 
zone  an  algebraic  grid  is  generated  with  grid  clustering  near  walls  atul  high  How  gradient 
regions.  The  computational  method  is  constructed  such  that  each  zone  is  considered  an  inde- 
pendc’iit  module  interacting  with  other  zones  bi'fore  or  after  the  informat  ioti  corresponding 
to  each  zone  is  updated  one  cycle.  Zonal  boundaries  are  transparent  to  the  (luwlield.  A 
typical  (i-zone  grid  used  for  computations  described  in  this  ix'port  is  designed  as  follows  (see 
f  igures  5  and  G):  grid  zone  1  covers  the  projectile  from  nose  to  base,  zone  '2  covers  tlx' 
aica  between  zone  1  and  the  inner  surface  of  the  sabot,  zone  1  covers  the  area  l)etween  the 
outer  surface  of  the  sabot  and  th<'  uppermost  extent  of  the  computational  domain,  zones  .a 
and  G  cover  the  projectile  and  sabot  ha.se  regions,  rr'spect i vely.  Zones  1  thin  G.  excluding 
zone  d.  extend  from  c)  =  0  to  60^'  in  tlx'  azimuthal  direction.  Cirid  zone  H  covc'is  the  ari'a 
between  the  sabot  and  the  azimuthal  extent  of  the  com|)ulat ional  domain.  I  lu'  (uitire  G 
zone  grid  consists  ol  HOU.OGO  nodes  aiul  ri'cjuires  10  million  words  of  memoiy  on  a  ('HA^  2 
su|)('rcomputer.  ('onverg<-d  si>lutions  re<|uir<' about  10  CPI'  hours. 
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3.  RESULTS 


I'iffuic  7  shows  tlu'  measured  (Selimidt  1981)  and  computed  pressure  disf ril)ut ions  over 
(he  projectile  and  sabot  surface  in  the  pitch  plane;  the  pitch  plane  (see  Figur<'  1)  bisects 
the  azimuthal  nlanform  of  the  sabot.  Three  sabot  components  are  modeled  with  sabot 
ba.ses  aligiie  with  the  projectile  base.  Ax/ D  =  0,  projectile  surface  and  sabot  inner  surface 
vertically  se|jarated  by  Ay / D  —  .7.5,  with  the  sabots  at  zero  angle-of-attack.  Laminar 
boundary  layer  modeling  was  einj>)(jyed,'  turbulent  solutions  are  described  elsewhere  (Nusca 
.•\ug.  1990.  Oct.  1990).  Computed  pressures  on  the  projectile  surface  agree  favorably 

with  the  magnitude  and  location  of  a  measured  pressure  peak  {x/D  ~  1.22)  as  well  as 
elevat('(l  pre.ssurcs  preceding  this  peak,  2  <  x/D  <  4.22.  The  location  of  this  pressure  peak 
corrc’sponds  to  the  termination  of  a  low  speed  flow  region  on  the  projectile.  Downstream  of 
the  iM'<'ss\ire  peak  the  agreement  between  computation  and  measurement  is  also  favorable. 
On  the  inner  surface  of  the  sabot,  numerical  .simulation  ad('quately  predicts  the  pressure 
h'vel  and  (rend  on  the  sabot  slant  surface,  2.75  <  x/D  <  3.91.  Prc'ssure  trends  on  the  rest 
of  the  sabot  s<*ction  agree  with  measurements  including  a  pressure  ris('  at  x /  !)  ~  5.5. 

.\usca  ( 1990,  1991 )  describes  further  results  obtained  for  the  wind  tunnel  model,  summa¬ 
rized  here.  For  cases  when  the  sabot  petals  are  close  to  the  projectile  {Ay/ D  <  .75),  a  low- 
speed  ( :\{  <  I )  recirculating  flow  pocket  forms  between  the  projectile  and  the  beveled  section 
of  the  .sabot  petals.  This  causes  a  strong  impinging  oblitiue  shock  on  the  projcctil(>  surface 
where  the  jrocket  forms  and  a  high  pressure  pulse  where  the  pocket  terminates.  ,A.s  the  sabot 
petals  di.scard,  a  normal  shock,  formed  at  the  leading  edge  of  the  sabot .  becomes  an  oblique 
shock  that  intersects  the  projectih*  surface  in  a  regular  reflection.  Inviscid  flow  simulations 
re<)uire  significantly  less  computer  time  by  excluding  the  viscous  terms  in  the  Navier-Stokes 
equations.  However,  the  inviscid  simulation  predicts  lower  pressures  on  the  projectile  and 
sabot  than  measured  or  predicted  by  laminar  and  turbulent  simulations.  Turbulent  calcu¬ 
lations  are  similar  to  laminar  for  the  low  Reynolds  number  wind  tunnel  data.  Comparison 
of  CFD  predictions  with  projectile  surface  data  measured  azimuthally  about  the  projectile 
agre<'  with  tlu-  tr'^nd  but  not  the  magnitude  of  these  pressures  (in  particular  the  pressure 
peak,  as  shown  in  Figure  7,  reduces  as  njeasured  azjmuthally  about  the  body).  Azimuthal 
grid  refinement  increases  the  level  of  agreement.  Computations  for  the  2D/axisymmetric 
eq\iivalent  of  three  sabot  petals  (i.e.  petals  joined  into  a  concentric  tube  with  the  i)rojectile 
cc'uterline)  are  computationally  inexpensive  but  result  in  flowfields  that  are  very  different 
from  the  .31)  case. 

Figures  S  through  17  show  comi)uted  laminar,  steady-state,  pressure  contours  in  the 
pitch  |)lane  for  the  forward  part  of  the  projectile/sabot  configuration  and  ten  stages  of  the 


programmed  discard  sequence.  Three  horizontal  lines  extending  from  i  j D  —  0  to  7.03  are 
zonal  grid  boundaries.  Large  flow  gradients  (e.g.  shock  waves)  are  indicated  by  clustering  of 
pressure  contour  lines.  Pressure  contour  levels  are  the  same  for  Figures  9-17,  1  <  P/Poo  ^ 
10,  AP/Poo  =  -5;  however,  for  Figure  8,  1  <  PI  Poo  <  ^0^,  API  Poo  —  1  due  to  higher 
stagnation  pressures. 

The  programmed  discard  sequence  shown  in  Figures  8-17  covers  four  vertical  displace¬ 
ments  of  the  sabot  inner  surface  with  respect  to  the  projectile  surface  {Ayl D)  and  six  sabot 
angles  of  attack  (with  respect  to  the  projectile).  The  projectile  was  assumed  to  be  at  zero 
yaw  with  respect  to  the  freeslream  and  the  Mach  number  was  constant  as  1.5.  Sinc«'  the  time 
during  which  the  sabot  petals  and  projectile  are  in  close  proximity  is  usually  short  (about  2 
nis  or  1.5  meters  from  the  gun),  the  assumption  of  constant  Mach  number  is  not  unreason¬ 
able.  This  quasi-steady,  programmed  simulation  ignores  the  flow  time  dynamics  and  does 
not  link  the  aerodynamic  forces  to  the  sabot  motion.  However,  such  a  simulation  serves  as 
a  prelude  to  computations  that  utilize  coupling  of  unsteady  aerodynamics  and  rigid-body 
motion. 

.'\s  seen  in  Figures  8-17,  the  sabot  generates  a  strong  series  of  shock  waves,  beginning 
as  a  detached  nearly-normal  shock  that  intersects  the  projectile  surface  as  a  strong  ol)iique 
shock,  and  ending  as  an  attached  oblique  shock  that  intersect  the  projc'Ctile  surface  in  a 
regular  reflection.  Flow  betw(>en  the  sabot  inner  surface  and  the  projectile  surface  Itegins  as 
a  choked  uearly-uniform  high  pressure  field  with  transition  into  reflected  shocks  (from  sabot 
back  to  projectile)  that  l)ecoiue  more  pronounced.  Beginning  with  Figure  12.  a  low  pressure 
recirculation  bubble  develops  on  the  sabot  inner  surface  extending  from  .i'll)  =  3.t)l  to  the 
u<-.\t  shock  irnping('ineut  on  the  sabot  surface.  Combinedwith  tlie  high  pressure  on  the  sabot 
be\eled  sectit.  (2.75  <  .v I D  <  3.9-1)  this  low  (pressure  region  j)rovides  a  force  coujde  that 
pi\^motes  .sabot  discard. 

lisiug  the  simulated  sabot  discard  .sequence  described  above,  the  corresponding  aerody¬ 
namic  forces  (lift  and  drag)  and  pitching  moment  can  be  computed.  This  is  accomplished 
by  integrating  the  sabot  surface  pressure  and  shear  stress  distributions  for  each  stage  of 
the  discard  .sequence.  The  sabot  mass  properties  are  used  to  compr.tc  vertical  and  hori¬ 
zontal  accelerations  which  are  assembled  in  a  table  as  functions  of  sabot  Ayl  I)  and  n.  A 
modified  point-mass  trajectory  model  is  used  to  compute  the  sabot  cc'uter  of  gravity  (CC) 
location  as  a  function  of  time  using  double-interpolation  from  values  in  the  tabh*.  l-'igure  18 
shows  a  comf)arison  between  the  sabot  CC  location  (both  in  the  axial  and  radial  directions) 
compute'd  using  the  AVCO  semi-enipirical  code  and  ‘he  present  simulation  using  CFD.  I'he 
pre.sent  pre-dictions  match  the  AVCO  values  for  early  time-s,  but  diverge  later  in  the  simu¬ 
lated  discard  event.  In  the  AVCO  simulation  sabot  discard  progresses  faster  than  predicted 


usiii^  the  current  met  hod.  I'iie  relatively  good  agre<*meiit  for  <'arly  times  in  the  discard 
event  may  he  a  resnll  of  the  sabot /(nojectile  interference  methods  incinried  in  th(“  A\  (’() 
corh'.  Keasons  lor  (liscrejrancics  in  the  predictions  at  later  tiriK's  are  still  under  investigation. 
One  possibility  is  that  the  Nevvtoniati  theory  used  to  predict  aerodynamic  forces  when  the 
sabot  is  not  in  clo.se  proximity  to  the  projectile,  re.sults  in  lift  and  drag  values  that  are  larger 
than  pre(licted  nsitig  CFI).  Iti  comparing  the  AVCO  i)rediction  to  that  using  Cl"  I),  several 
points  should  be  noted.  Both  methods  used  the  same  sabot  gt'ometry  atid  tnass  properties, 
freest  ream  How  conrlitions  and  assumed  a  symmetric  discard,  l^otli  methotls  are  (piasi-str'ady 
in  nature,  usitig  a  database  of  stearly  aerodynamic  force  predictions  to  simulatr'  a  dynamic 
ev('nt.  Howi'ver,  the  source  of  the  acModynamic  data  is  very  rlifhnent  between  the  codes  (sec 
■‘Introduction"  for  a  discussion  of  tlie  AVCO  code).  By  virtue  of  the  rapid  aerodynamic 
methods  incorporated  into  the  .AV'CO  corle,  a  much  larger  aerodynamic  force  and  moment 
database  is  available.  The  trajectory  time-integration  stej)  for  th(“  AVCO  cod*'  was  much 
smaller  that  that  used  in  the  prr'sent  study. 

I'igitre  1!)  shows  the  piojectile/.sabot  configuration  of  tli<>  Army  .M8G')  anti-iatik  roiitid. 
riu'  coidiguratioti  has  beeti  altered  .somewhat  in  order  to  simplify  cotnpntat  ional  grid  gener- 
atioti.  1  hese  alterations  are  also  illustraterl  in  Figure  It).  Fix'  sabot  was  located  .7.')  calibers 
above  the  |)roj('Ctile  (1  caliber  =  d8mm)  and  at  zero  atigh'-of-attack.  I'ignre  20  shows  the 
com|)nfatiotial  grid.  A  simulated  sabot  discard  sccjnence  like  that  n.sed  for  th('  wind  tunnel 
model  is  in  progress.  Figure  21  shows  the  laminar  flow  pressure  contours  for  the  M86.o.  The 
Ucynolds  nnml)er  for  this  flow  is  6.6  million  per  meter  and  the  freestrearn  Mach  number  is 
4..n. 


4.  CONCLUSIONS 

CFD  solutions  of  the  .fl)  Navier-.Stokes  equations  have  been  applied  to  the  aerodynamics 
of  symmetric  sabot  discard.  ,\  steady  simulated  sabot  discard  sequence  using  fi.xed  sabot 
locations  (with  respect  to  the  projectile)  reveals  shock/shock  aiid  shock/boundary-layer  in- 
tf'ractions  in  tl)e  flowfield.  Tlie  fieestream  Mach  number  was  1.")  and  laminar  boundary  layer 
modeling  was  employed  for  Reynolds  number  6.6  million  per  meter.  NunH'rical  simulations 
have  also  been  performed  using  Reynolds  number  89  million  per  meter  and  flows  with  tnr- 
bidencr'  modeling  (Nusca  Aug.  1990).  ddie  steady-state  approach  that  uses  |  redetermined 
sabot  |>osition.s  has  lead  to  enhanced  understanding  of  the  discard  event,  serving  as  a  prelude 
to  computations  that  utilize  coupling  of  unsteady  aeio  ly  .;amics  and  rigid-body  motion.  .A 
technirpre  for  the  integration  of  surface  pressures  and  shear  stre.s.s  was  developed  for  the  wind 
tunnel  model  sabot. 
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Numerical  mesh  geiieratioii  lor  the  solution  of  complex  flowfields  about  realistic  projec¬ 
tile/sabot  configurations  may  be  greatly  simplified  by  the  use  of  unstructured  (i.e.  finite- 
element  like)  grids.  Figure  12  shows  the  planar  view  (i.e.  slice  through  the  j)itch  plane  of 
tlie  projectile/sabot)  of  a  typical  unstructured  grid  for  the  Army  M829  sabot.  Solution  of 
the  Filler  equations  on  unstructured  grids  is  being  accomjilished  by  Ciiakravartliy  (1991). 
Work  on  unstructured  grids  and  moving  grid  zones  will  eventually  lead  to  a  more  realistic 
simulation  of  the  discard  event. 
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Figure  1. 


Figure  2 


Photograph  of  typical  kinetic  energy  Jong-rod  projectile  in 
during  three-petal  sabot  discard. 


.  Shadowgraph  of  typical  kinetic  energy  long-rod  projectile 
during  four-petal  sabot  discard. 
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Figure  3.  Schematic  of  wind  tunnel  model  in  the  pitch  plane  ((f)  =  0. 180°). 
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Figure  7.  Laminar  flow  pressure  distributions  for  projectile  and  sabot  surfaces 
in  the  pitch  plane  (<i>  =  0, 180°),  A.r/Z)  =  0.  Ay/D  =  .75.  a  =  0°. 
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Figure  11.  Laminar  flow  pre.s.surf'  rontours  in  tlic  pitch  i)lan<’  (o  =  0.  ISO'  )  fur 
1 1)  =  0,  l\]il D  =  0.75.  o  =  0'’. 


17 


Projectile 


Figure  12.  Laminar  How  pressure  contours  in  the  pitch  plane  (o  =  0.  180'^)  foi 
Aj'//;  =  0.  Av/Z)  =  0.75,  a  =  2°.  ~ 


■>  T: 


Projectile 


Figure  13.  Laminar  flow  pressure  contours  in  tlie  pilch  plane  {o  =  0,  180°)  for 


Ax /I)  =  0.  Ay/D  =  0.75.  o 
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Figure  14.  I.aininar  flow  pressure  contours  in  the  pitch  plan<’  {o  —  0.  18( 
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Figure  15.  Laniinar  flow  pressure  contours  in  the  pilch  plane  ( o  =  0.  180  ')  fm 
Ax /I)  =  0.  Ai/fl)  =  0.75,  o  =  8®. 
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Figure  16.  Laminar  flow  pressure  contoius  in  the  pitch  plane  (o  =  0.  ISO'")  fui 
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Figure  17.  Laminar  flow  pressure  contours  in  the  pitch  iilane  (o  =  0.  ISO")  for 


Ax/ D  =  0,  Ay/D  =  1.0,  o 
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Figure  19.  M865  projcctile/sabot  configuration.  Solid  line  is  actual  geometry. 
Dashed  line  is  computational  geometry. 
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Figure  21.  Laminar  flow  pressure  contours  in  the  pitch  plane  ((/>  =  0, 180°) 
for  M865  sabot,  Ai/D  =  .957,  Ay//)  =  .75.  a  =  0°. 


Figure  22.  Unstructured  grid  (pitch  plane  view)  for  M8‘29  projectilc/sabot. 
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LIST  OF  SYMBOLS 


cal 

C’p 

Cv 

I) 

t 

i\a,n 

M 

Pr 

P 

'ft 

T 

t 

V 

u.  i\  w 

II' 

■'  -i/, 

A.r 

A/y 


=  caliber,  D 

=  specific  heat  capacity,  constant  p 
=  specific  heat  capacity,  constant  volume 
=  diameter  of  projectile  cylinder  section 
=  specific  total  internal  energy 
=  flux  vectors  (Eq.  1) 

=  Mach  number 
=  Prandtl  number 
=  pressure 
=  heat  transfer  rate 
=  specific  gas  constant 
=  temperature 
=  time 

=  mean  streamwise  velocity 
=  cartesian  velocity  components 
=  dependent  variable  vector  (Eq.  1) 

=  cartesian  coordinates 

=  horizontal  distance  between  projectile  and  sabot  ba.ses 
=  vertical  distance  between  projectile  and  sabot  surfaces 


Clreek  Symbols 


o  =  sabot  angle  of  attack  wrt  the  projectile 

O  =  ratio  of  specific  heats,  Cp/c„ 

C  =  transformed  coordinate 

r/  =  transformed  coordinate 

p  =  molecular  viscosity 

i  =  transformed  coordinate 
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!>  ~  (Iciisily 

^jr-  ^:iy  =  iioiinal  st icss  U-nsors 

T  —  t lansforniocl  time 

’■j-,/,  Tjz,  Ty,  =  sliear  stress  tensors 

(f>  —  aziimitiial  angle,  0  and  180“  for  pitch 

plane,  clockwise  looking  downstream 


Subscripts 


/ 

.r,/y,  c 

c 

oc 


=  turbulence  quantity 
=  denotes  spacial  components 
=  (,'-direct ion  trausforiti  coefficient 
=  //-direction  transform  coefficient 
=  i^-direction  transform  coefficient 
=  freestream  quantity 
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ATTN:  F.E.  Wolf 
BaralKK).  W1  53913 

3  Olin  Ordnance 

ATTN:  E.J.  Kirschke 

A.F.  Gon/alcz 

D. W.  Worthington 
P.O.  Box  222 

St.  Marks.  FL  32355-0222 

1  Olin  Ordnance 

ATTN:  H.A.  McElroy 
lOIOl  9ih  Street.  North 
,St.  Petersburg.  FL  .3.3716 

1  Paul  Gough  AsscK'iatcs,  Inc. 

ATTN:  P  S.  Gough 
1048  South  St. 

Portsmouth,  NH  0.3801-542.3 
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No.  of 
Copies 


No.  ()1 

C'of'ies  ( )i  Li:iiii/;ilioii 

1  i'liysics  lnicrnalion;il  Library 
A  TIN:  H.  Wayne  Wampler 
F*.0.  Box  .sOlO 

San  Leaiulro.  CA  94577-05^ 

2  Prineelon  Combuslion  Research 

Laboralories,  Ine. 

ATTN:  N.  Mer 

N.A,  Messina 
Prmeeion  Coiporaie  Pla/a 
1 1  Deerpark  Dr..  Bldg  IV.  Suite  119 
Monmouth  Junction.  NJ  08K.s2 

Rockwell  International 
Roeketdyne  Division 
.ATTN:  BA08. 

J.  Flanagan 
J.  Gray 
R.B.  Edelman 
66.^.^  Canoga  Avenue 
Canoga  Park.  CA  91.70.V27()3 

2  RiKkwell  International  Science  Center 
ATrN:  Dr.  S.  Chakravarthy 

Dr.  S.  Piilaniswamy 
1049  Camino  Dos  Rios 
P.O.  Box  lOH.s 
Thousand  Oaks.  C A  91.160 

1  Science  Applications  International  Corp. 

ATTN;  M.  Palmer 
2109  Air  Park  Rd. 

Albuquerque.  NM  87106 

I  Southwest  Research  Institute 

ATTN;  J.P.  Ricgcl 
6220  Culcbra  Road 
P.O.  Drawer  28.M0 
San  Antonio.  TX  78228-0.*)  10 

1  Sverdrup  Technology,  Inc. 

ATTN;  Dr.  John  Deur 
2001  Aerospace  Parkway 
Brook  Park.  OH  44142 

3  Thiokol  Corporation 
Elkton  Division 
ATTN;  R.  Wilier 

R.  Biddle 
Tech  Library 
P.O.  Box  241 
Elkton.  MD  21921-0241 


Organi/ation 

1  Veritay  Technology.  Inc. 

ATTN;  E.  Fisher 
484.5  Millcrspori  Hwy. 

East  Amherst.  NY  14.501-0305 

1  Universal  Propulsion  Company 

ATTN;  H.J.  MeSpadden 
25401  North  Central  Avc. 

Phwnix.  AZ  85027-7837 

I  SRI  International 

Propulsion  Sciences  Division 
ATTN;  Tech  Library 
333  Ravenw(HHl  Avenue 
Menlo  Park.  CA  94025-3493 

Aberdeen  Provmi;  Ground 

I  Cdr.  USACSTA 

ATTN;  STECS-PO/R.  Hendricksen 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-204 _ Date  of  Report  Septeaber  1993 

2.  Dale  Report  Received _ 

3.  Docs  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so.  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  Slate,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


Department  of  the  Army 


OFRCIAL  BUSINESS 


i 

1 


BUSINESS  REPLY  !VL\IL 
niisT  cKss  Bam  it  mi.  w 

Postage  *tU  be  p»>d  by  addiessee 


Director 

U.S.  Army  Research  Laboratory 
ATTN;  AMSRL-OP-Cl-B  {Tech  Lib) 
Aberdeen  Proving  Ground,  MD  21005-5066 


